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Abstract 

We compute the energy levels of a 2D Hydrogen atom when a constant 
magnetic field is applied. With the help of a mixed-basis variational method 
and a generalization of virial theorem, which consists in scaling the wave func- 
tion, we calculate the binding energies of the IS, 2P~ and 3D~ levels. We 
compare the computed energy spectra with those obtained via a generalization 
of the mesh point technique as well as the shifted 1/N method. We show that 
the variational solutions present a very good behavior in the weak and strong 
magnetic field regimes. 

In the past few years a great effort has been devoted to the study of two-dimensional 
atoms. The technological advances in semiconductor processing has enabled to manu- 
facture nanostructures where the electrons are confined to two dimensional (quantum 
wells), one dimensional (quantum wires) and zero dimensional (quantum dots) mo- 
tions (HI, ?, |, [3|, |, U and references therein). Because of the reduced dimensionality, 
new and interesting transport and optical properties can be observed. In order to 
explain the experimental data, a large number of theoretical investigations have been 
performed H [7|, ||. 

Probably the best known quantum well configuration consists of regions of GaAs 
which act as wells of conduction electrons, separated by regions of Gax^AL-As which 
act as barriers. The application of magnetic field perpendicular to the quantum wells 
is expected to provide further band-structure data and binding energy information. 
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The importance of the theoretical calculations demands reliable analytical and nu- 
merical techniques. In this study we are motivated to analyze the binding energies of 
the ground and excited states of a hydrogenic donor in a 2D electron gas. Here, with 
the help of a scaled mixed variational method, we compute the energy levels of a 2D 
Hydrogenic atom in the presence of a magnetic field of arbitrary strength,we compare 
our results with those obtained by numerical computation and with the shifted 1/N 
method 
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The Hamiltonian describing the electromagnetic interaction between a conduction 
electron and a donor impurity center when a constant magnetic field B is applied 
perpendicular to the x-y plane can be written as 
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where we have written the vector potential A =y(— y, x, 0) = ^fe v in the symmetric 
gauge. The coupling constant 7, which measures the ratio between the magnetic en- 
ergy and Coulomb energy, is defined as 7 = e 2 h 3 B/(ce 3 m* 2 ) where m* is the effective 
mass, e the dielectric constant of the host material, V 2 is the two-dimensional Lapla- 
cian, and L z is the angular momentum operator —ihd/dcf) with eigenvalue hm.The 
units of energy are given in terms of the effective Rydberg constant TZq = m*e 4: /2h 2 e 2 
and the effective Bohr radius a* = h 2 e/m*e 2 , respectively. 
The substitution 



*(r) 

reduces the Schrodinger equation H^f 
differential equation 



Efy to the following second order ordinary 
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Exact solutions of eq. (j3|) cannot be expressed in closed form in terms of special 
functions. There are analytical expressions for the energy for particular values of 7 
and m [|l|, 13, 14 . Besides numerical and perturbation methods [15, different 
techniques have been used in order to obtain the eigenvalues E in equation (Rl), in 



particular the two-point Pade approximation [17, 18 



Recently, using the shifted 1/N expansion, Mustafa [[LI] , [19| has computed the IS, 
2P~ and 3D~ energy levels for a 2D donor impurity in the presence of a magnetic 
field. The shifted 1/N and Pade methods are very powerful in computing the energy 
spectra but fail in giving a reasonable simple wave function and general criteria for 
the lower bound of energy values. In the present article we discuss the problem using 
a scaled two-terms mixed-basis variational approach. We compare our results with 



those obtained using the Schwartz p(J interpolation technique, as well as the shifted 
1/N method; showing that our results fit very well to those computed numerically 
for any range of values of the magnetic strength B. 
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In order to apply the variational method to our problem ||21| , we look for a trial 
wave function. Since equation @ reduces to the Hydrogen atom equation when 
7 = 0, we can consider as a basis for 7 << 1 the Hydrogen wave functions Since 
< ^ h \H\ >> E, we obtain a lower bound of the energy for small values of the 
parameter 7. The solution of equation ([|) when 7 = is 

<Mp) = An^e-^/^+HyH^ 2 H J P ) (4) 

(1/2 + n p + \m\) 



where D m n is a normalization constant, and L(a,b,x) are the Laguerre polynomials ||23||. 
Consequently the energy spectrum in the zero-field limit takes the form 

1 

Eh = ~ (1/2 + n p + \m\y (5) 

Conversely, for large values of 7, a good trial basis is that of the spherical oscillator. 
In this case the solution of equation (H) has the form 



Vos C (r) = C<^e- w2/ V m| L(n p , |m| , |p 2 ) (6) 

and in the high- field limit, the energy levels are 

E sc = l(2n p + \m\+m + l). (7) 

If we attempt to apply the variational method using the hydrogen atom basis, 
we will obtain good agreement with the accurate results for small values of 7, but 
this approach fails for large 7 even if we consider a many term basis. An analogous 
situation occurs when we use the oscillator basis, in which case we obtain a good 



agreement for large 7 but the convergence is very slow for small values of 7 |22| . In 
order to solve this problem, we propose a mixed-basis approach. The idea is to use 
as trial function, for any quantum level, a linear combination of the form 

^ = C H ^ H + Co^Osc (8) 

where and ^osc are the corresponding hydrogen and oscillator wave functions 
associated with the quantum level in question; Co and ch are constants to be calcu- 
lated. It is worth noticing that our mixed basis is not orthogonal under the inner 
product J °° ^>i^> jpdp. We proceed to minimize the expectation value \H\ with 
the normalization condition, | \&) = 1 

Applying the variational approach to the basis coefficients q, we reduce our prob- 
lem to that of solving the matrix equation 

[{% \H I - A {% I Cj = 0, (9) 
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whence after substituting the Hamiltonian (JT|) into (^), and taking into account the 
separated differential equation (|3]) we obtain the matrix equation, 



E H + \A + mi -\ (E 0sc -X)S-2C 
(Eqsc — A) S — 2C Eq — 2B — A 



Cl 
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where 
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where the lowest value of A will be the energy of the level. Since we obtain a sec- 
ond order equation for A, we can analytically compute the energy eigenvalues and 
eigenvectors. The advantage of this approach is twofold. First, because of variational 
approximation, the eigenvalues satisfy the inequality A > E and therefore we have a 
lower bound for our energy levels. Second, we obtain a relatively simple expression 
for the normalized eigenf unctions. 

Instead of restricting ourselves to the the energy values obtained from ([)]) via the 
mixed-basis variational method, we can improve the accuracy of the energy with the 
help of the virial theorem. Using this technique, after re-scaling the radial parameter 
r — > £r we have that 

9{p) - *(£p) 



and the matrix terms associated with p n and the second derivative become 

(*< \p n \ *;> - r n i(£p) n i *i(&0) = r n (*< \ P n \ 

d 



d 



dp 2 
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(12) 
(13) 
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Then, with the help of the virial theorem for a Coulomb and a radial oscillator 
potential, 

Eh - 
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-(*h\T h \*h) = -{*h\V h \*h) 

Eosc = 2 (^Osc \Tosc\ ^Osc) = 2 ($Osc \Vosc\ ^Osc) 

where T is the kinetic energy term, the re-scaled Hamiltonian reduces to 
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with 

J = / (— h m ^H^oscjpdp. (17) 

Jo ctp ap 

In order to compute the quantum energy values, we proceed to compute the lowest 
eigenvalue of the matrix equation given by 

( i2n(0-A H 12 (£)-\S 
{ H 12 (0-\S fl 22 (£)-A 





which takes the form 



, _ flu + H 22 — 2SH\ 2 , . 



(flu + fl 22 - 2Sfl 12 ) 2 - 4(1 - S 2 ){H U H 22 - H 
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2(1 -S 2 ) 

Since the components of fl depend on the scale parameter £, A is a function of £, 
and the optimal value for £ can be obtained computing the minimum for A(£), which 
reduces to solving the algebraic equation 

dX , , 

ST - (20) 

Then, we can substitute the value of £, obtained after solving (|20|) , into ( |TI|) and 
obtain an analytic expression for the scaled wave functions. As illustration of the 
mixed-basis method, we present two different expansions of the ground state of th 
2D Hydrogen Hamiltonian ([]]). We plot the energy against 7' = 7/(7 + 1) as the 
horizontal scale. 

For comparison, the numerical computations of the energy spectra associated with 
equation (|3|) are carried out with the help of the Schwartz method [20| which is a 



generalization of the mesh point technique for numerical approximation of functions. 
This method gives highly accurate results given a thoughtful choice of the reference 
function. For equation (|[) we chose as the interpolation function 

/(P) = E/™ , , (21) 

where 

u(p) = sin[7r(p//i) 1/2 ]. (22) 

Here r m is a zero of u(p), a m is a zero of its derivative, and h is the step of the 
quadratic mesh. Here, as illustration of the mixed-basis method, we present two 
different expansions of the ground state of the 2D Hydrogen Hamiltonian 

Figures 1, 3, and 5 compare the variational results with those obtained numerically 
and with the help of the shifted 1/N method. 
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Figure 1: Energy of the IS state as a function of 7'. The thin solid line is obtained by 
numerical methods; the dotted line is obtained using the mixed IS hydrogen and IS 
oscillator bases variational method. The thin broken line is obtained via the scaled 
variational method with the IS hydrogen and oscillator bases. The thick broken line 
is obtained with the help of the shifted 1/N method 
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Figure 2: The figure shows the difference between the numerical result for the IS en- 
ergy spectrum and the energy values computed with the help of the mixed variational 
(thin broken line), scaled variational (thin solid line), and the shifted 1/N method 
(thick broken line) 
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Figure 3: Energy of the 2P~ state as a function of 7'. The thin solid line is obtained 
by numerical methods; the dotted line is obtained using the mixed 2P~ hydrogen 
and 2P~ oscillator bases variational method. The thin broken line is obtained via 
the scaled variational method with the 2P~ hydrogen and oscillator bases. The thick 
broken line is obtained with the help of the shifted 1/N method 
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Figure 4: The figure shows the difference between the numerical result for the 2P~ 
energy spectrum and the energy values computed with the help of the mixed vari- 
ational (thin broken line), scaled variational (thin solid line), and the shifted 1/N 
method (thick broken line) 
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Figure 5: Energy of the 3D~ state as a function of 7'. The thin solid line is obtained 
by numerical methods; the dotted line is obtained using the mixed 3D~ hydrogen 
and 3-D" oscillator bases variational method. The thin broken line is obtained via 
the scaled variational method with the 3D~ hydrogen and oscillator bases. The thick 
broken line is obtained with the help of the shifted 1/N method 
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Figure 6: The figure shows the difference between the numeric result for the 3D~ en- 
ergy spectrum and the energy values computed with the help of the mixed variational 
(thin broken line), scaled variational (thin solid line), and the shifted 1/N method 
(thick broken line) 
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It is easy to see that even for a two term mixed basis a good fitting is obtained in 
the weak and strong field regimes. One of the bases gives a reasonable good fitting in 
the intermediate region. A better fit is obtained with the help of the scaled variational 
method. We also have that for the 2P~ and 3D~ states the mixed-basis variational 
approach gives very good results. 

Examining Figs. 1, 3 and 5, it is not obvious which technique gives the most 
accurate results for computing the hydrogen energy levels. However, Figures 2, 4 and 
6 show that the shifted 1/N method always gives results below the the numerical 
energy levels. Among the three analytic methods, figures 2, 4 and 6 show that the 
scaled variational method gives the most close to numerical results even for large 
values of 7. 

It would be interesting to apply the mixed-basis technique as well as the scaled 
variational method for the 2D Hydrogen problem when relativistic effects are not 
negligible. This will be the object of a forthcoming publication. 
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